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! ABSTRACT 

We present the first cluster catalogue extracted from the UKIRT Infrared Deep Sky 



00 

o 



Survey Early Data Release. The catalogue is created using UKIDSS Ultra Deep Survey 
infrared J and K data combined with 3.6 ^m and 4.5 /im Spitzer bands and optical 
BVRi'z' imaging from the Subaru Telescope over 0.5 square degrees in the Subaru 
| XMM-Newton Deep Field. We have created a new cluster-detection algorithm, based 

"--»« . on the Friends-Of- Friends and Voronoi Tessellation methods, which utilises probability 

' distribution functions derived from a photometric redshift analysis. We employ mock 

catalogues to understand the selection effects and contamination associated with the 
Q ■ algorithm. The cluster catalogue contains 13 clusters at redshifts 0.61 ^ z ^ 1.39 with 

?-) | luminosities 10L* ^ L tot ^ 50L*, corresponding to masses 5 x 10 13 M & M c i uster & 

3 x 10 14 M for = 75h. The measured sky surface density of - 10 deg~ 2 for 

high-redshift [z = 0.5 — 1.5), massive (> 10 14 M ) clusters is precisely in line with 
theoretical predictions presented by Kneissl et al. (2001). 

Key words: galaxies: clusters: general - catalogues - methods: analytical - methods: 
data analysis - surveys - galaxies: photometry - galaxies: high-redshift 



1 INTRODUCTION 

Distant galaxy clusters are used in a wide range of cosmolog- 
ical and astrophysical contexts. For instance, cluster studies 
illuminate how dark matter haloes collapse and large-scale 
structure forms and evolves. Comparisons of the evolution of 
the cluster mass-function to that predicted by iV-body sim- 
ulations (Evrard et al. 2002) or the Press-Schechter formal- 
ism (Press & Schechter 1974) and its variants place strong 
constraints on cosmological parameters such as the mass 
density (Qm) and the amplitude of the mass fluctuations 
in the early Universe (as) (e.g. Eke et al. 1998). Further, 
high-redshift clusters can help constrain feedback processes 
caused by star-formation and Active Galactic Nuclei (e.g. 
Silk & Rees 1998). Unfortunately few z > 1 clusters have 



been identified, and the vast majority of these are from X- 
ray surveys (e.g. using XMM-Newton; Stanford et al. 2006). 

X-ray selected cluster surveys (e.g. Mullis et al. 2005) 
have succeeded in finding clusters with luminous intra- 
cluster media. Some follow-up observations exploiting the 
Sunyaev-Zel'dovich (SZ) effect have been successful (e.g. 
Jones et al. 2005) and comprehensive SZ surveys are un- 
derway (Kneissl et al. 2001). Optical selection of clusters is 
complementary to these studies (Gilbank et al. 2004) as it 
is based on other cluster properties, e.g. optical vs. X-ray 
luminosity or the SZ decrement of the cluster gas. An added 
advantage of using optical data is the possibility of employ- 
ing a photometric redshift analysis to derive the distance to 
the clusters and gain a three-dimensional perspective. Opti- 
cal cluster searches however have been stymied at high red- 
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shifts by the fact that the 4000A break (ubiquitous for the 
early- type, red galaxies dominant in clusters) shifts out of 
the /-band into the near-infrared. Recent developments have 
seen the advent of large near-infrared cameras like the Wide 
Field Infrared Camera (WFCAM) on the United Kingdom 
Infrared Telescope (UKIRT). WFCAM is now undertaking 
the UKIRT Infrared Deep Sky Survey (UKIDSS, Lawrence 
et al. 2006). 

In this letter we use the UKIDSS Ultra Deep Survey 
(UDS) Early Data Release (EDR, Dye et al. 2006) to find 
clusters in the redshift range 0.5 ^ z ^ 1.5. The clus- 
ters are found by applying adaptations of two cluster selec- 
tion methods: the Voronoi Tessellation technique (Ebeling 
& Wiedenmann 1993) and the Friends-Of-Friends method 
(Huchra & Geller 1982). We describe the data and the de- 
termination of photometric redshifts in §2, and the cluster- 
detection algorithm and simulations using mock catalogues 
in §3. The cluster catalogue is presented in §4 and in §5 we 
summarize our conclusions. Throughout this letter we as- 
sume Ho = 71 km s~ 1 Mpc~ 1 , fi M = 0.27, and Q A = 0.73. 
All magnitudes quoted are in the AB-system. 



2 DATA AND PHOTOMETRIC REDSHIFTS 

We use three sources of data: near-infrared J and K data 
from the UDS EDR (Foucaud et al. 2006); 3.6/xm and 4.5/um 
bands from the Spitzer Wide-area InfraRed Extragalactic 
survey (SWIRE, Lonsdale et al. 2005); and optical BVRi'z' 
Subaru data over the Subaru XMM-Newton Deep Field 
(SXDF, Furusawa et al. in prep.). The area of our cluster 
survey is 0.5 sq. deg. with 34.2° < RA < 34.8° and -5.4° < 
Dec. < -4.6.° (02:16:48.00 < RA < 02:19:13.20, -05:25:12.0 
< Dec < —04:38:17.1). The galaxy catalogue includes the 
objects with a detection in i' , J and K and to exclude stars 
we impose a criterion of SExtractor stellarity index < 0.8 
in i' and K (e.g. Bertin & Arnouts 1996). The depth of the 
catalogue is limited by the UDS EDR 5a magnitude limits of 
J^AB,iim = JAB.iim = 22.5. We derive photometric redshifts 
by fitting a spectral energy distribution (SED) template to 
each object's photometry using the Hyperz code (Bolzonella 
et al 2000). The galaxy templates are generated with the 
stellar population synthesis code GALAXEV (Bruzual & 
Chariot 2003) and cover a range of different star formation 
rates with timescales r from 0.1 to 30 Gyr. We adopt a flat 
prior for galaxy luminosity up to a maximum of L — 10L* 
(assuming a passively evolving elliptical galaxy) in the ob- 
served A'-band and calculate the marginalised posterior red- 
shift probability functions. Our photometric redshift cata- 
logue comprises 19300 objects in the range 0.1 ^ z ^ 2.0. 

2.1 Reliability of the photometric redshifts 

To test the reliability of the photometric redshifts we run 
Hyperz with the same templates on spectroscopically iden- 
tified objects in the SXDF (Yamada et al. 2005; Simp- 
son et al. 2006). The spectroscopic and photometric red- 
shifts are consistent, with a photometric redshift error of 
<t z = 0.08 over 0.5 5j z ^ 1.5 and no obvious systematic off- 
set. However, an examination of the number of galaxies in 
the photometric-redshift catalogue as a function of redshift, 
N(z), reveals spikes in the redshift distribution. The shape 



of the distribution is constant across the field which indi- 
cates the spikes are not caused by large-scale structure but 
by aliasing effects inherent to photometric redshift meth- 
ods employing many more templates than data points (e.g. 
Rowan- Robinson 2003). A serious concern is the detection of 
spurious clusters due to this focussing in redshift. We there- 
fore create a second redshift catalogue using four empirically 
derived SEDs from Coleman, Wu & Weedman (1980), sup- 
plied with the Hyperz code. The new galaxy redshift proba- 
bility functions are broader; the spikes disappear yielding a 
smoother N(z) distribution. However, a comparison of the 
spectroscopic redshifts reveal that the redshift is systemati- 
cally underestimated at z <; 0.5, with an offset of Az = 0.15 
at z ~ 1. 

In summary, an examination of the accuracy of the pho- 
tometric redshifts reveals a sensitivity to the templates used 
in the photometric redshift estimation. Redshifts quoted in 
this letter use the Bruzual & Chariot templates as the spec- 
troscopic sample show these to be most accurate. However 
to avoid false cluster detections due to the spikiness in the 
N(z) distribution we limit ourselves in this letter to the clus- 
ters that are isolated in both photometric redshift catalogues 
and are therefore considered to be robust detections. 



3 THE ALGORITHM 

The cluster-detection algorithm is described in detail in van 
Breukelen et al. (in prep.). Here we present only a brief 
outline of the procedures we employ. 

Two common problems of optically selected cluster 
samples are the detection of spurious clusters, and the 
chance projection of fore- and background galaxies into the 
real clusters. This is caused by selection biases inherent to 
any detection algorithm and the fact that photometric red- 
shift probability functions (z-PDFs) can take any form, ex- 
hibiting errors that may be many times larger than the red- 
shift range of the cluster. To reduce these difficulties we (i) 
use two substantially different cluster-detection methods to 
minimize the number of false detections, and (ii) utilize the 
full z-PDFs rather than the best redshift estimate with an 
associated error. For each method we sample the z-PDFs to 
create 500 Monte-Carlo (MC) realisations of the 3-D galaxy 
distribution. These are divided in slices of Az = 0.05 1 in 
which the cluster candidates are identified. Since the er- 
ror on the photometric redshift is larger than the width of 
the redshift slices, each cluster-candidate is typically found 
in several adjoining slices. We determine the final cluster- 
redshift by taking the average of the redshifts at which the 
cluster occurs, weighted by the corresponding number of 
MC-realisations. We assign a 'reliability factor' F to each 
of the clusters which is the total fraction of MC-realisations 
in which the cluster is detected. The algorithm uses the 
Voronoi Tessellation technique (VT), and the Friends-Of- 
Friends method (FOF) to detect cluster candidates in the 
redshift slices. One of the principal advantages of the VT 
method is that it is relatively unbiased as it does not look 
for a particular source geometry (e.g. Ramella et al. 2001). 

1 When the analysis was re-run with the alternative photometric 
redshift catalogue (Sec 2.1) this was changed to Az = 0.1 to 
account for the broader z-PDFs. 
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The field of galaxies is divided into Voronoi Cells, each con- 
taining one object: the nucleus. The reciprocal of the area of 
the VT cells translates to a local density. Overdense regions 
in the plane are found by fitting a function (see Kiang 1966) 
to the density distribution of all VT cells in the field; clus- 
ter candidates are the groups of cells of a significantly higher 
density than the mean background density. The Friends-Of- 
Friends (FOF) algorithm groups galaxies with a smaller sep- 
aration than a projected linking distance Dn n k ('friends'). 
When spectroscopic data is available, the 'friends' are also 
subject to a linking velocity AV ^ Mink (e.g. Tucker et al., 
2002). With photometric redshifts only galaxies within the 
same redshift slice are linked (see Botzler et al. 2004). If a 
group comprises a number of galaxies larger than n m i n it is 
considered a cluster candidate. 

We create mock catalogues that mimic the real data to 
quantify the various selection effects and expected contam- 
inants. These simulations also lead to estimates of the total 
stellar luminosity of the clusters found in the real data. 

3.1 Simulated Cluster Catalogues 

We create catalogues with a galaxy distribution randomly 
placed 2 in the field with 0.1 ^ z ^ 2.0. The galaxy lumi- 
nosities and number densities are determined by the A"-band 
luminosity function of Cole et al. (2001), with the simpli- 
fying assumption of passive evolution with formation red- 
shift Zform = 10. A detection limit of K < 22.5 is imposed 
to match the 5-cr limit of the UDS EDR data. The num- 
ber of galaxies as a function of magnitude in each mock 
catalogue is entirely consistent with the number counts in 
the UDS catalogue to this limit. We superimpose simulated 
clusters of total mass M = 0.5,1.0,2.0,10,20 x 10 14 M Q , 
assuming M J^ e = 75/i (Rines et al. 2001) and applying a 
cluster if— band luminosity function (Lin et al. 2004). The 
galaxies are spatially distributed within the cluster accord- 
ing to an NFW profile (Navarro, Frenk & White, 1997) with 
a cut-off radius of 1 Mpc. For the simulated clusters of mass 
> 1.0 x 10 15 Mq this is smaller than the virial radius; how- 
ever we do not find such high mass clusters in our data (see 
Section 4) and the effect of the cut-off is therefore negligi- 
ble. The clusters are placed at redshifts z = 0.1, 0.2, 2.0. 
Each combination of mass and redshift is represented in ten 
randomly created catalogues. We offset the redshift of all 
galaxies randomly with a factor equal to the errors in the 
real data. Each galaxy receives a z-PDF taken from the data 
by selecting an object with corresponding redshift. 

The VT and FOF methods each have two free param- 
eters. For FOF these are the linking distance in proper co- 
ordinates, -Dunk, and the minimum number of galaxies in a 
cluster, n m i n . Guided by Botzler et al. (2004) we experiment 
with values between 0.125 Mpc ^ -Diink ^ 0.175 Mpc, and 
3 ^ n m i n 5. For VT the parameters are the maximum 
probability of an overdensity being a background fluctua- 
tion, pbg, and the lower limit on the cell density, f m - m . We 
follow the method of Ebeling & Wiedenmann (1993) and set 
Pbg to 10%. For / m i n we try values of 1.2 — 2.2, where / = 1.0 
is the mean cell density of the field. We use the parameters 
that optimize the algorithm's performance: Dimk = 0.175 

2 We neglect clustering of both the background and the clusters. 



Mpc, n m m = 5, and / m i n = 1.74. Since both methods use 
different measures to isolate clusters (galaxy density in VT 
versus separation in FOF) the false detections in both do 
not coincide. Therefore by cross-correlating the output of 
the two methods the spurious sources due to biases in the 
algorithms disappear, which reduces the contamination to 
chance galaxy groupings. There is no obvious bias due to 
cluster morphology (for a full treatment of biases and con- 
tamination see van Breukelen et al in prep.). 

We compare the recovered cluster galaxies to the in- 
put galaxies of the mock clusters. VT tends to include all 
galaxies in a large area around the cluster core (see for exam- 
ple the left-hand panel of Fig. . The number of recovered 
cluster members, N sa \, in any cluster is sensitive to the lo- 
cal field density. By contrast, the galaxy members recovered 
by FOF are more centrally concentrated; the total number 
per cluster is consistent throughout the random realisations 
of the catalogues. Thus we use both methods to detect the 
clusters and only FOF to calculate iVg a i. This is done by 
taking all galaxies that occur in the cluster at > 15% of the 
MC-realisations in which the cluster itself is detected. The 
galaxies that appear in a smaller fraction of MC-realisations 
are very likely to be interlopers from different redshifts. Cal- 
culating iV ga i for all cluster-masses at all redshifts yields 
functions of N gB ,\ vs. z for constant mass or total luminos- 
ity. At z > 1.5 only clusters of masses M ^ 1.0 x 10 15 M Q 
or L > 170L* can be detected with ~ 50% completeness. 
At redshifts of z < 0.5 our field of view is ^ 10 times the 
typical cluster size of 1 Mpc. Hence we limit our redshift 
range to 0.5 ^ z ^ 1.5, allowing the detection of clusters of 
luminosity L i; 10L* at z ~ 0.5 to L <; 170L* at z ~ 1.5. 

4 RESULTS: THE CLUSTER CATALOGUE 

Our final cluster catalogue comprises all clusters found by 
both VT and FOF with a reliability factor of F > 0.2 to ex- 
clude false detections. We calculate the redshift by averaging 
the cluster redshift given by the VT and FOF methods. This 
results in 14 clusters at 0.61 ^ z ^ 1.39; TableQlists all the 
clusters with their positions and properties (full details of 
the cluster galaxies are available on request). The redshift 
error given in Table combines the errors given by VT and 
FOF and reflects the range of redshift slices in which the 
cluster is detected. Although the cluster-finding method is 
sensitive to clusters down to a luminosity of ~ 10L* we 
are limited by the relatively large systematic errors from 
the photometric redshifts. The 13 clusters 'above the line' 
in Table Q were also detected in our alternative photomet- 
ric redshift catalogue (see Section 2.1), although sometimes 
only at the F{ { and F vt > 0.1 level; thus these are judged 
robust detections. Only one cluster ('below the line' in Ta- 
ble was not recovered with the alternative photometric 
redshift catalogue. Fig. shows three cluster examples. On 
the left a _K"-band image, a combined Bz'K image in the 
middle, and on the right a colour-magnitude (CM) diagram. 
For illustration a modelled red sequence is overplotted on the 
CM-diagram (using GALAXEV elliptical templates of Zf orm 
= 10, t = 0.5 Gyr and a slope derived from Kodama & Ari- 
moto 1997). A clear red sequence can be seen for the top two 
clusters; the third is unclear because of the small number of 
cluster galaxies, although four of the seven galaxies detected 
in both VT and FOF lie near the predicted red sequence. 
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Figure 1. Clusters 1, 6, and 13 spanning the redshift range of our detections. Left: _ft"-band images; the large circle shows a 1 Mpc region 
around the cluster; the blue squares and red circles are cluster members as given by VT and FOF respectively. ID 1: the green cross 
denotes an extended X-ray source. ID 6: the green arrows point out galaxies with z Bpec = 0.87. Middle: Bz' K images of the central 1 
Mpc region. Right: colour-magnitude plots of the clusters: the colour is the reddest filter shortward of 4000A (restframe) minus K. The 
crosses are all galaxies within the central 1 Mpc region, otherwise the symbols are the same as in the left-hand panel. The grey band is 
the modelled red sequence. 



Three clusters (no. 1, 2, 11) consist of two concentra- 
tions of galaxies separated by ^ 1 Mpc in the FOF method, 
but as one cluster in VT: these could be merging clusters. We 
determine JV ga i with K < 22.5 (corresponding to the com- 
pleteness limit) in the same way as for our simulated clusters 
(see Section 3.1); this allows us to derive an approximate 
total luminosity to the cluster by interpolating between the 
lines of constant total luminosity in the iVgai — Z plane found 



in our simulations. We find our clusters span the range of 
10L* L tot ^ 50L*; assuming -jfj^ = 75h (Rines et al. 
2001) this yields 0.5 x 10 14 M <, M clustor <, 3 x 10 14 M Q . 
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Table 1. The cluster catalogue ordered by redshift. The position (J2000) is the centroid of the cluster galaxies. -/V ga i is the number of 
cluster galaxies detected by FOF with K < 22.5; Ltot is the estimated total luminosity of the cluster and mass the inferred total cluster 
mass; F( a { and F v t are the reliability factors given respectively by FOF and VT. The possible spectroscopic redshift is given as £spec ; 
clusters marked with a star [*] are shown in Fig. Q All clusters above the line are robust; ID 14 was not recovered with the alternative 
photometric redshift catalogue, t Extended X-ray-detected cluster (Kolokotronis et al. 2006); z spC c from Geach et al. (in prep). J Two 
galaxies in each of the clusters are found at this redshift (Yamada et al. 2005). ft Four galaxies in the cluster are found at this redshift 
(Yamada et al. 2005). Potential cluster galaxies have z spec =0.874, 0.96, and 1.095; super-position effects may be exaggerating the 
richness of this cluster, ttf Redshift of a QSO possibly associated with the cluster (Sharp et al. 2002). 



5 CONCLUDING REMARKS 

Our cluster catalogue comprises 13 clusters at redshifts 0.61 
^ 2 ^ 1.39 with estimated luminosities 10L* Si Ltot Si 50L*, 
corresponding to masses 5 x 10 13 M© M c i ustcr ^ 3 x 
10 Mq . Considering just the clusters with total mass esti- 
mates significantly > 10 14 M (6, 10, 11, 12, 13), this rep- 
resents a sky surface density of ~ 10 deg -2 for high-redshift 
(z = 0.5 — 1.5) clusters. This is in quantitative agreement 
with the predictions of low-density (Qm ~ 0.3) cosmolog- 
ical models, e.g. those presented by Kneissl et al. (2001). 
Comparing with their Figs. 5 & 8, we see that, within the 
obvious limitations of small-number statistics, the real clus- 
ters in the SXDF have the same abundance, redshift and 
estimated mass distributions as are predicted by theory. 

Clearly spectroscopic observations are essential to con- 
firm the reality of the clusters. From spectroscopic redshifts 
in the literature (see Table 0, we have tentative confirma- 
tion of the photometric redshifts for 6 out of 9 clusters at 
2 > 0.75. In the near future near-infrared multi-object spec- 
trometers on 8-metre class telescopes will provide opportu- 
nities for spectroscopic follow-up of high-redshift clusters. 

As we note only one of our cluster candidates is as- 
sociated with a published extended X-ray-detected cluster 
(Table 0, we plan to exploit existing deep XMM-Newton 
data (Ueda et al. in prep.) and upcoming SZ surveys to de- 
termine the gas contents of these clusters. 
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